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ABSTRACT 

We present the results of two (101 ks total) Chandra observations of the 
z = 0.778 optically selected lensing cluster RCS022434-0002.5, along with weak 
lensing and dynamical analyses of this object. An X-ray spectrum extracted 
within R2500 (362 kpc) results in an integrated cluster temperature of 5.1^0 5 
keV. The surface brightness profile of RCS022434-0002.5 indicates the presence 
of a slight excess of emission in the core. A hardness ratio image of this object 
reveals that this central emission is primarily produced by soft X-rays. Further 
investigation yields a cluster cooling time of 3.3 x 10 9 years, which is less than half 
of the age of the universe at this redshift given the current ACDM cosmology. A 
weak lensing analysis is performed using HST images, and our weak lensing mass 
estimate is found to be in good agreement with the X-ray determined mass of 
the cluster. Spectroscopic analysis reveals that RCS022434-0002.5 has a velocity 
dispersion of 900 ± 180 km s" 1 , consistent with its X-ray temperature. The core 
gas mass fraction of RCS022434-0002.5 is, however, found to be three times lower 
than expected universal values. The radial distribution of X-ray point sources 
within R 2 oo of this cluster peaks at ~ O.7R200, possibly indicating that the cluster 
potential is influencing AGN activity at that radius. Correlations between X-ray 
and radio (VLA) point source positions are also examined. 

Subject headings: cosmology:observations — X-rays:galaxies:clusters — galaxies:clusters:general 
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Introduction 



Clusters of galaxies which exhibit grav itationa l lensing have been actively sou ght af- 
ter since their discovery in the mid-1980s (jSoucaill Il988l ; iLynds &: Petrosianl Il988l ). This 
interest is due to the numerous scientific opportunities that these objects provide. Be- 
cause they magnify sources which lie at very high redshifts, they allow us to investigate 



form i ng galaxies in the early epochs of the universe in great deta il (e.g., lEbbels et al. 



19961 ; iFranx et al. 1 119971 ; iPettini et al. I 1200(1 ISwinbank et al. 1 120031 ). Gravitational lens 
ing can also be used as an independent probe of the mass distribution of cluste rs (e.g., 



Kneib et al. 1993 ; Tyson. Kochanski. fc delF Antonio 1998 ; Broadhurst et aL~ll2005l ). In ad- 
dition, clusters with multiple lensed sources conta in information on the underlying cosmol- 
ogy of the universe (IGolse. Kneib. fc Soucailll200ll ). Until recently the cumulative sample of 
strong l ensing clusters was relatively small and mainly restricted to objects at redshifts of 
z < 0.5 JWu et al. Ill998l : I Williams. Navarro, fc Bartlemannlll999f ). 

The Red-Sequence Cluster Survey (RCS; Gladders fc Yeel 200ol, 2005t) has recently dis- 
covered eight strong lensing clusters at high redshift ( Gladders et al. 20031 ). One of the most 

striki ng objects discovered via the RCS survey is RCS022434-0002.5 at z = 0.778 (IGladders. Yee. fc Ellingsc 



20021 ). This cluster possesses a rich array of lensing arcs that originate from several back- 
ground sources, one o f which lies at a (spect roscopically confirmed) redshift of z = 4.88, 
anoth er at z = 3 . 66 (ISwinbank et al. 1 120071 ) , and a radial arc was found to be at z = 



1.050 (ISand et al. 112005] 



RCS022434-0002.5 (hereafter RCS0224-0002) was observed twice with the Chandra X- 
ray Observatory, initially for 15 kiloseconds, and more recently for 89 kiloseconds, the deepest 
single X-ray follow-up observation of an RCS cluster to date. Here we present our analysis 
of these observations. From the data we derive a spatially resolved surface brightness profile, 
cluster temperature, and mass. In addition, we compare high resolution X-ray flux and hard- 
ness images of the cluster core with observations in the optical and radio. Weak lensing and 
dynamical analyses are performed and compared to X-ray results. A point source analysis 
of this field is also presented. 



Unless otherwise noted, this paper assumes a cosmology of Ho = 70 km s _1 Mpc -1 , 
Qyi = 0.3, and fi^ = 0.7. Using this cosmology, 1" = 7.42 h^Q kpc at the cluster redshift. In 
addition, all error bars represent 68% confidence levels. 
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2. Chandra Observations and Data Reduction 

RCS0224-0002 was observed by the Chandra X-ray Observatory on 15 November 2002 
for 12,327 seconds, and 9 December 2004 for 89,202 seconds. The aimpoints of the obser- 
vations were located on the S3 chip of the Advanced CCD Imaging Spectrometer (ACIS) 
CCD array. Both observations were conducted in VFAINT mode, and reprocessed accord- 
ingly. Upon obtaining the data, aspect solutions were examined for irregularities and none 
were detected. Charged particle background contamination was reduced by excising time 
intervals during which the background count rate exceeded the average background rate by 
more than 20%. Bad pixels were removed from the resulting event files, and they were also 
filtered on standard grades. After cleaning, the available exposure times were 12,051 seconds 
and 88,973 seconds, providing a total exposure of 101,024 seconds. 

Combined images, instrument maps and exposure maps were created in the 0.29-7.0 
keV band using the CIAO 3.3.0.1 tools DMCOPY and MERGE_ALL, binning the data by 
a factor of 2. Data with energies above 7.0 keV and below 0.29 keV were excluded due 
to background contamination and uncertainties in the ACIS calibration, respectively. Data 
between 0.29 keV and 0.6 keV were included despite questionable low energy calibration for 
the purpose of increasing the overall signal-to-noise ratio of the observation. 

Flux images were created by dividing the image by the exposure map, and initial point 
source detection was performed by running the tools WTRANSFORM and WRECON on the 
0.29-7.0 keV image. These tools identify point sources using "Mexican Hat" wavelet functions 
of varying scales. An adaptively smoothed flux image was created with CSMOOTH. This 
smoothed image was used to determine the X-ray centroid of the cluster, which lies at an RA, 
Dec of 02:24:34.161, -00:02:26.44 (J2000). The uncertainty associated with this position can 
be approximated by the smoothing scale applied to this central emission (2.5"). Determined 
in this manner, the position of the X-ray centroid is 4.5 ± 2.5" (33.4 ± 18.6 kpc) from 
the brightest cluster galaxy (BCG) of RCS0224-0002, falling within the ran ge of typical 



offsets between BCGs and X-ray isophotal centers in X-ray selected samples (Patel et al. 



20061 ). Figure [T] shows an HST image of RCS0224-0002 with adaptively smoothed X-ray flux 



contours overlayed. 



3. Surface Brightness 

A radial surface brightness profile was computed from a combined, point source removed 
image and exposure map over the range 0.29-7.0 keV in 1" annular bins, then multiplied by 
a factor of (1 + z) 3 to correct for cosmological dimming. We find an excess of emission within 
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the central ~ 4" of the cluster. This central excess does not appear to be an AGN, since 
it is extended (~ 4x wider than the PSF) and no hard point source was detected at that 
location. To test the validity of this central excess, a second surface brightness profile was 
computed in 2" annular bins. 

Both profiles were then fit with both single and double (3 models. A single (5 model 
takes the form: 



r 2\ -3/3+! 



I(r) = I B + I 1 + , 1 



where Ib is a constant representing the surface brightness contribution of the background, 
Io is the normalization and r c is the core radius. 

A double j3 model is represented as: 



where each component has fit parameters (J n , r n , j3 n ). 

The central excess does persist in the more highly binned surface brightness profile, yet 
the reduced x 2 s of the single vs. double f3 model fits are almost identical in both cases, and 
there is no clear statistical preference between the two models. (3 model fits to both the 
1" and 2" bin surface brightness profiles are shown in Figure [2] and details of all four fits are 
given in Table [TJ Unless otherwise noted, the single f3 model fit to 2" bin data is used in all 
subsequent analysis. 



4. Spectral Analysis and R 



2500 



Using the centroid position indicated in Section (2J a point source removed spectrum 
was extracted from each observation in a circu lar region wit h a 300 h^ 1 kpc radius. These 
spectra were then jointly analyzed in XSPEC (|Arnaudlll996l ). using weighted response ma- 
trices (RMF) and effective area files (ARF) generated with the CIAO tool SPECEXTRACT 
and the Chandra calibration database version 3.2.2. Backgrounds were extracted from the 
aimpoint chip in regions free of cluster emission. The resulting spectra, one extracted from 
the short observation and one from the long observation, containing 183 and 1651 source 
counts respectively. Both spectra were grouped to include at least 20 counts per bin. 
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The spectra were jointly fit with a single temperature spectral model including fore- 

groun d absorption. The absorbing column was fixed at its measured value of 2.91 x 10 20 cm -2 (IDickey & Loci 



19901 ) . and the metal abundance was fixed at 0.3 solar (jEdge &: Stewartlll99ll ). Data with 
energies below 0.29 keV and above 7.0 keV were excluded from the fit. 

The results of these fits, combined with the single f3 model parameters from Section [3j 
were then used to estimate the value of R2500 f° r this cluster. This is accomplished by 
combining the equation for total gravitating mass 



M (< ) kT(r)r f 5 In p ^ 5 In T\ 
G/j>m p \5 In r 5 In r J ' 

(where pm p is the mean mass per particle) with the definition of mass overdensity 



M tot (r A ) = -irp c (z)r 3 A A, (4) 

where z is the cluster redshift, and A is the factor by which the density within r& exceeds 
Pc(z), the critical density at z. Here p c (z) is given by p c (z) = 3HqE(z) 2 /8nG, where E(z) = 
[Cl m (l + z) 3 + ^a] 1//2 . These equations are then combined with the density profile implied 
from the f3 model (assuming hydrostatic equilibrium, spherical symmetry, and isothermality) 



-3/3/2 



Pgas{r) = Po 



1 + 



(5) 



resulting in the equation 



(Ettori 


2000; 


Ettori et al. 


2004) 



After the initial estimation of -R25oo ; additional spectra were extracted from within 
that radius, and spectral fitting was performed again. This iterative process was continued 
until the fitted temperature and calculated value of i?25oo was consistent with the extraction 
radius. The final estimate of -R2500 is 3621^ h^ 1 kpc (48.8"), with an extrapolation to i? 2 oo 
of 1.47[ti4i hfQ 1 Mpc. The best fit temperature of the spectra extracted from within this 
region was T x = 5.1^5 keV, with a reduced x 2= 0-75 for 77 degrees of freedom. Using this 
spectrum, the unabsorbed bolometric luminosity of RCS0224-0002 within -R2500 is L x = 



3pkT 
GW(4/3)vrp c (z)r c 2 A 



(6) 
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2. 3^0 3 x 10 44 erg s 1 . Figure [3] is a plot of these spectra with the best fitting model overlaid 
on the data. 



5. Hardness Map 

Images, instrument maps and exposure maps were created for both a soft X-ray band 
(0.5-2.0 keV) and a hard band (2.0-8.0 keV) by the method described in Section [2j These 
files were then used to create flux images in the two bands. 

Using the smoothing scale produced in the analysis of the entire energy band, smoothed 
flux images were created in each of the two bands. The hard band image was subtracted 
from the soft band image, and the result of that calculation was divided by the smoothed flux 
image of the entire energy band (i.e., Is I ~| H )■ Contours of the resulting hardness ratio image 
are overlayed on a VLA 1.4 GHz radio image in Figure HI Soft emission is concentrated in a 
region containing both the X-ray centroid and the BCG, a radio galaxy to the southwest of 
the X-ray centroid (Figure H]). There is no clear emission peak that coincides with the BCG, 
a radio galaxy at z — 0.7773. 

6. Mass Determinations 

6.1. Gas Density and Cooling Time 

The gas mass density distribution of a cluster of galaxies whose surface brightness profile 
is well described by a (3 model can be shown to follow Equation This relationship is 
dependent upon the assumptions of hydrostatic equilibrium, spherical symmetry, and gas 
isothermality. The central density in this equation can be determined by an expression 
relating the observable cluster X-ray luminosity to gas density. This calculation requires the 
experimental determination of both gas temperature and emission measure, which can be 
performed in XSPEC. 

With the results of spectral analysis (Section HJ) and using the best fitting surface 
brightness parameters for the inner luminosity peak (Section [3]), a central gas density of 
7.6 ± 0.09 x lO" 3 cm" 3 is obtained for RCS0224-0002. 

Recall that the characteristic time that it takes a plasma to cool isochorically through 
an increment of temperature 5T can be written 
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5t 



cool 



3 k 
2 n A(T) 



5T, 



(7) 



( ISarazmlll988l ). where n is the electron density, A(T) is the total emissivity of the plasma 
(the cooling function), and k is Boltzmann's constant. In the case of isobaric cooling, 3/2 is 
replaced by 5/2. 

Using this formula combined with the cooling function appropriate to a 0.3 solar abun- 
dance plasma, the cooling time of RCS0224-0002 is 3.3 x 10 9 y ears. The age of the univers e 
in our adopted cosmology at z — 0.778 is 6.8 x 10 9 years (e.g., iThomas fc Kantowskil 12004 ). 



Therefore the cooling time of the innermost region of RCS0224-0002 is less than the 
age of the universe at z = 0.778. This, combined with the luminosity excess in the central 
few arcseconds, as well as the existence of a soft X-ray peak at the same location, indicates 
that this cluster may contain a proto-cooling core. The observation does not contain enough 
source photons to obtain spatially resolved temperature information for this cluster, therefore 
longer exposures will be required to completely characterize the central emission. 

At least one other study has d e tected a high redshift cluster with the signatures of 
a small cooling core. iGrego et al. I (120041 ) in their analysis of a Chandra observation of 
MS 1137. 5+6625 (at z = 0.783) detect a small central surface brightness excess, and an 
excess of soft emission in the core of this cluster. They also calculate a central cooling time 
that is shorter than the age of the universe at that epoch. 



6.2. X-ray Mass Estimate 

Using the results of the (3 model fits (Section EJ) and spectral fit (Section Hj), along with 
Equation [5] and the equation of hydrostatic equilibrium (Equation [3]), it can be shown that 



M, n ,(< r) 



3(3 kTr {r/r c Y 



G jj,m p l + (r/r c ) 



2 ' 



Using this equation (with \i = 0.62), gas mass and total mass were calculated out to R2500 
(~ 362 /ifo 1 kpc) and R 200 (~ 1-47 Mpc) for RCS0224-0002. One sigma errors on central 
density, gas mass, total mass, and gas mass fraction were determined through statistically an- 
alyzing a random sampling (N = 10, 000) of values drawn from the 68% confidence intervals 
of relevant parameters, assuming that their errors are uncorrelated. Mass determinations 
resulted in M 25 oo, t ot = 1.6 ± 0.2 x 10 14 h^ 1 M , M 25 oo,gas = 0.05 ± 0.005 x 10 14 h TO 5/2 M , 
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M 200 ,tot = 8.3 ± 0.8 x 10 



ii 



h^n Mm and M 



SQO^as = 0.24 ± 0.02 x 10 14 h 70 5/2 M . The mass- 



1.2 ±0.05 x 10 14 hfo 1 



temperature relation of IVikhlinin et al.l (120061 ) predicts M 2 5oo,tot 
for z = 0.778 and Tx = 5.1 keV, which is within 2cr (25%) of our observed value. 



M 



© 



The core gas mass fraction of RCS0224-0002 is /, 



g,2500 



0.031 ±0.005 h 



-3/2 
70 • 



This value 



is significantly lower than that d etermined via the s ame method for a sample of 14 medium 



redshift X-ray selected clusters ( IHicks et al. 1 120061 ) and is three to four times lower than 



both theoretically and observationally expected "universal" values (jEvrardlll997l ; lAllen et al. 



2004J). This is a phenomenon that is also seen in groups, poor clusters, and other high red- 



shift clusters ( IDeir Antonio. Geller. &: Fabricantl Il995l: ISanderson et al. I l2003t ISadat et al. 



20051 ). a s well as high reds hift cluster simulations (jEttori et al. 1120041 ) and additional RCS 



clusters (IHicks et al. 1 120071 . in prep). 



6.3. Weak Lensing Analysis 



The central region of RCS0224-0002 was observed in the FQ0QW and F8UW filter 
with the WFPC2 camera onboard the Hubble Space Telescope, with the center of the cluster 
located on WFC3. The observations consist of 2 orbits in each filter, yielding total integration 
times of 6600s in FQ06W and F8UW. 

The exposures in each passband were split into two sets (by the orbit the data were 
taken), resulting in two images. These two images were desi gned to be offset by 5.5 



pixels, to allow the construction of an interlaced image (e.g., Ivan Dokkum et al. 1 11999 
Hoekstra. Franx. fc Kuijkenl 120001 ) , which has a sampling that is a factor \/2 better than 
the or iginal WFPC2 images. This same approach was used by iHoekstra. Franx. fc Kuijken 
(120001 ) in their study of MS 1054-03 and we refer to this paper for more details. For the anal- 
ysis presented here we omitted data from the Planetary Camera because of the the brighter 
isophotal limit. 



The next step is to use the interlaced images for our weak lensing ana 



ber of modifications described in Hoekstra et al. 



ysis. 



The analysis 

technique is based on the one devel oped bvlKaiser. Squires, fc Broa dhurst ( 1995) with a num- 



1998 



(119981 ) and IHoekstra. Franx. fc Kuiiken 



(2000 ). We note that this im plementation has been tested extensively (e.g. 



Hoekstra et al. 



Heymans et al. 1120061 ) and has been shown t o be a ccurate at the few percent level. We 



follow the procedure outlined in IHoekstra et al. I (119981 . 2000) which results in a catalog of 
ellipticities for the faint galaxies that we use in our lensing analysis. These shapes have been 
corrected for PSF anisotropy and the size of the PSF. 



The measurement of the shape of an individual galaxy provides only a noisy esti- 
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mate of the weak gravitational lensing signal. We therefore average the measurements of 
a large number of galaxies as a fu nction of distance to the cluster center. As discussed in 
Hoekstra. Franx. fc Kuijkenl (120001 ) we weight each object with the inverse square of the un- 
certainty in the distortion, which includes the contributions of the intrinsic ellipticities and 
the shot noise in the shape measurement. 

Figure [5] shows the resulting average tangential distortion as a function of distance from 
the cluster center, which was taken to be the location of the brightest cluster galaxy. We 
selected 462 galaxies with 22.5 < F8UW < 25.5 and FQOQW - F8UW < 1.5 (i.e., bluer 
than the cluster early types) as source galaxies. We detect a significant lensing signal. The 
lower panel of Figure [5] is the signal, gx, when the phase of the distortion is increased by 
7r/2. If the signal observed in the upper panel is due to gravitational lensing, g x should 
vanish, as is observed. The solid line is the best fit singular isothermal sphere model which 
yields an Einstein radius of 4.7 ± 0.9 arcseconds. 

To relate the observed lensing signal to an estimate of the mass requ i res an estimate of 
the mean source redshift distribution. As in lHoekstra. Franx. fc Kuijkenl (120001 ). we use the 



results from lFernandez-Soto. Lanzetta. fc Yahill (119991 ) which are based on the Hubble Deep 
Fields. This allows us to determine an equivalent line of sight velocity dispersion using the 
best fit isothermal sphere model. We obtain a value of a = 777l 7 g km s _1 . 

Alternatively we can fit an NFW model to the measurements. We assume that the 
density profile i s only a functio n of th e virial mass and use the relation between concentration 
and mass from iBullock et al. I ( 20011 ). This yields a virial mass of 5. 41 4 / \ x 10 14 M and a 
corresponding M 250 o = 1-3^^6 x 10 14 M (and r 2 5oo = 375 kpc). This result is in excellent 
agreement with both the dynamical and X-ray analysis results. 



As shown in iHoekstra et al. I (120021 ) mass estimates based on single HST pointings tend 
to be biased low, because of substructure in the cluster center. However, the distribution 
of galaxies in RCS0224-002 is relatively regular and does not show significant substructure. 
Furthermore, because of its high redshift, the lensing signal can be measured sufficiently far 
out to be less affected by substructure. 



6.4. Dynamical Analysis 

Optical multi-object spectroscopic observations of galaxies in this field were performed 
by several telescopes and instruments; dates and observational parameters are listed in Ta- 
ble |2j Both galaxies and gravitational arc images were targeted, over a field of view of 5 
arcminutes. Most cluster galaxies had magnitudes of R~ 20-22, though a few galaxies as 
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faint as R~ 24 were identified via their [Oil] emissio n lines. Several slits targeting cluster 



arcs were tilted or curved (see iGladders et al. 1 120031 ). but these were not used for cluster 



galaxy redshifts. Overall, 39 galaxy redshifts at 0.70 < z < 0.85 were measured by the 
wavelength of individual lines or via cross-correlation against galaxy spectra, with a range 
of spectral types. We find no significant systematic differences in the velocities measured 
by the different telescopes and instruments, either in velocity centroid, or in the 3 objects 
whose redshifts were measured accurately by more than one telescope. 

Figure [6] shows a histogram of galaxy velocities. There are two secondary peaks close 
to the cluster, one at z~ 0.750 and another at z~ 0.805. Both of these comprise galaxies 
spread over the field of view and are not associated with an easily identified isolated clump. 
There are substantial gaps between them and the cluster peak at z=0.778, and hence we do 
not identify these galaxies as cluster members. 

Correcting in quadrature for an estimated uncertainty of 100 km s _1 for individual 
velocity measurements, we find that the cluster is located at a mean redshift of 0.7777 and 
has a rest-frame velocity dispersion of 900 ± 180 km s _1 , based on 24 clust er member s. This 



estimate is very close to that expected from a standard a-T x relationship (jVoit et al. 112003 



Xue fc Wull2000l ). where an X-ray temperature of 5.1 keV predicts a velocity dispersion of 



872 km s _1 . Our velocity dispersion measurement is also consistent with b o th we ak lensing 
analysis (16.31) and the value of 920 km s _1 reported by ISwinbank et al. I (120071 ) based on 
their strong lensing analysis of this cluster. 



Point Sources 



7.1. X-ray Source Population 



Studies of the AGN populat ion in clusters of gal axies are an important probe of galaxy 
evolution in these environments (jMartini et al. 1120021 ) . While AGN can be difficult to detect 
at optical wavelengths due to obscuration, they are much more easily seen in the X-ray, 
where they comprise the dominant contribution to the hard X-ray point source population. 
In this section we perform an X-ray point source analysis of the Chandra observation of 
RCS0224-0002. 

Using the hard and soft band images and exposure maps referred to in Section [5j point 
sources in these bands were identified using the CIAO tools WTRANSFORM and WRECON. 
Only sources with a detection significance above 3a were retained for further analysis. In 
the case of faint source s, our numbers may represent a lower limit, due to both off-axis 
degradation of the PSF (IKim et al.ll2004j ) . and the possible masking of very faint sources by 
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diffuse X-ray emission in the cluster center. 



Using PIMMS (Portable Interactive Multi-Mission Simulator), assuming an absorbed 
powerlaw spectrum with a hydrogen column density fixed at the galactic value and a pow- 
erlaw index of 1.6, X-ray fluxes were calculated for each source after correcting for exposure 
variations across the chips. LogiV-logS' information was then calculated for each energy 
band, inclusive of all source s within R?n n . Th ese data were then compared to the logiV- 



\ogS backg round values o f iMoretti et al.l (120031 ). The results of these comparisons with la 
error bars (IGehrelsl 119861 ) are shown in Figure [7J The data are generally consistent with 
the measured X-ray background, however the background data do not extend to very low 
flux sources, and there may be a very slight excess of low-flux soft X-ray sources within 
R-200- "Cosmic variance", however, can become significant at low-flux as well, showing 
variations in source coun ts as high as 3.9a for 2 .0 — 8.0 keV sources with fluxes below 



< 1 x 10 15 erg cm 2 s 1 ( IBrandt fc Hasingerl 120051 ) 



A list of point sources within R200 is given in Tables[3]and TableHJfor the 0.5-2.0 keV and 
2.0-8.0 keV bands, respectively. These tables include both flux and luminosity information 
(also achieved with PIMMS, assuming a redshift of z = 0.778) for each source. 

To compare the dis t ributi on of sources in the RCS0224-0002 field to those reported 
in iRuderman fc Ebelingi (120051 . based on a statistical study of ~ 50 clusters), radial his- 
tograms containing the number of sources per square megaparsec (assuming a redshift of 
z = 0.778) were created from data in each energy band (Figure [H]). These plots indicate 
a si gnificant excess of point sou rces within R200, particularly in the 0.6 — 0.8 R200 bins. 
The IRuderman fc Ebelingi (120051 ) study was conducted in the 0.5-2.0 keV energy band, and 
our soft band data are roughly consistent with their results. 



7.2. Radio Comparisons 

RCS0224-0002 was observed with the National Radio Astronomy Observator y's (NRAO) 
Very Large Array (VLA)Q at 1.4 GHz for a total of 15 hours over three days (jWebb et al.l 



20051 ). Using this data, an investigation of the coincidence of X-ray point sources with 1.4 



GHz sources was performed. Figure [9] consists of the 1.4 GHz radio image of RCS0224-0002 
with 0.2-5.0 and 2.0-8.0 X-ray point source regions overlaid. Obvious correlations are noted 
in Tables |3] and IU and can be seen in Figure M None of the detected X-ray point sources 



1 The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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[in either energy band) is coincident with the cluster BCG. 



RCS0224-0002 was also observed at both 



Common-User Bolometer Array ([Holland et al 



:50jjm and 850/zm with the Submillimeter 
19981. SCUBA) in 2001 and 2002. The re- 



sults of these observations are reported in IWebb et al.l (120051 ). and have resulted in five 
submillimeter source detections in the RCS0224-0002 field. Only one of the five sources is 
coincident with an X-ray detected point source, smm-RCS0224.2 (source 6 in Table [3] and 
source 3 in Table HJ). 



8. Discussion and Summary 

This paper describes the results of the most in depth X-ray analysis of an optically 
selected high redshift RCS cluster to date, the 101 ks of Chandra observations of the lensing 
cluster RCS0224-0002 at z = 0.778. Imaging analysis reveals that the peak of X-ray emission 
is slightly displaced (4.5 ± 2.5") from the center of the l ensing arcs. This d isplacement falls 



well within the common range of X-ray centroid offsets (jPatel et al. 1 120061 ) . 



Single temperature spectral fitting of the X-ray emission within a 362 kpc (R2500) 
radius indicates that this cluster has an ambient temperature of 5.llo5 keV. The surface 
brightness profile of RCS0224-0002 exhibits a small central peak of emission which can be 
modeled via the inclusion of a second (3 component. The f3 model used in our analysis has a 
core radius of 213^}3 kpc and (3 value of 0. 88^0^4 . Using the results of both temperature 
and surface brightness fitting, R 2 oo is determined to be 1.47i$; 07 h^Q Mpc for this cluster. 

Hardness ratio maps were created from soft (0.5-2.0 keV) and hard band (2.0-8.0 keV) 
flux images. A peak in soft X-ray emission is seen to coincide with the central peak in surface 
brightness. This peak is elongated toward a 1.4 GHz VLA radio source which is centered on 
the BCG of this cluster, as seen in the HST image (Figure HJ). 

Using the results of both spectral and surface brightness fitting, we have calculated a cen- 
tral density for RCS0224-0002 of 7.6 x 10~ 3 ±9 x 10~ 4 cm -3 , which results in a central cooling 
time of 3.3 x 10 9 years. Without spatially resolved temperature information we are unable to 
confirm the presence of a cool component in the center of RCS0224-0002, however the cen- 
tral surface brightness excess, corresponding soft peak of emission, and lack of a central hard 
point source support the idea that RCS0224-0002 may contain a small proto -cooling core. 



Simil ar evidence has been detected in another z ~ 0.8 cluster, MS1137.5+6625 (jGrego et al. 



2004J). These two clusters are the highest redshift systems in which central cooling has yet 



been detected (albeit indire ctly) , and their exis tence provides evidence that massive clusters 



are in place at early times (jMullis et al. 112005 ) 
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Mass estimates of RCS0224-0002 result in total masses of M 



2500.tot 



10 



14 h 7 Q M and a gas mass of M 2 5oo,gas 



0.050 ± 0.005 x 10 14 h 



-5/2 



TO 



©■ 



= 1.6 ± 0.2 x 
This estimate 



of total gravitating mass is consistent within 68% errors with weak lensing mass estimates 
of M 250 o = x M©. Similarly, we find a rest-frame velocity dispersion for this 

cluster of 900 ± 180 km s" 1 , which is cons istent with that predicted by the a — Tx rela- 



tionship (IVoit et al. 1120031 ; Kue fc Wull2000f), as well as bo th weak (16. 3p and strong lensing 
estimates of the velocity dispersion (ISwinbank et al. 1 120071 ). The consistency of these mass 
estimators is an indication that RCS0224-0002 is a relatively relaxed system. 

An investigation of the point source population in the field of RCS0224-0002 indicates 
a possible slight excess of low-flux soft energy sources within R 20 o of the cluster center. The 
high energy source population of this cluster is generally consistent with the measured X-ray 
background. Both measurements of point source density represent lower limits because of 
a selection effect which makes it more difficult to detect point sources close to the cluster 
center (due to the surface brightness of the cluster). Radial distributions of these sources 
exhibit a significant peak within R 2 oo in both the hard and soft bands, possibly indicating 
that the cluster potential is influencing AGN activity at that radius. Spectral information 
on cluster member galaxies is needed, however, to test t hat hypothesis. The r a dial p rofile of 
soft X-ray point sources is consistent with the results of iRuderman fc Ebelina (120051 ). X-ray 
point sources within the field of RCS0224-0002 have been compared to 1.4 GHz VLA flux 
information, and coincident sources are listed in Tables [3] and HI 



The core gas mass fraction of this cluste r, .f g ,25oo 



0.031±0.005/i 



-3/2 
70 ' 



is more than three 

times lower than expected universal values (jEvrardlll997l ). T his is a phenomenon that is also 



1995 




Sanderson et al. 


2003 


tions ( 


Ettori et al. 


200^ 


[). It 



20051 ) . as well as high redshift cluster simula- 



mass of this object was converted into stars (jVikhlinin et all2006l ; lNagai. Kravtsov fc Vikhlinin 
20071 ); however additional analysis is required to confirm or rule out this possibility. In ad- 
dition, the AGN component of RCS0224-0002 may provide a significant contribution to the 
entropy of this system, however it is likely that other non-gravitational processes (i.e. galaxy 
formation) are also involved. 
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Fig. 1. — RCS0224-0002 HST Overlay, logarithmic X-ray contours are overlaid on an 
6600s HST image of RCS0224-0002 taken with the F606W filter. X-ray countours were 
created from an adaptively smoothed Chandra 0.29-7.0 keV flux image, and have values 
between 1.25 and 0.3 xl0~ 8 counts sec -1 cm -2 arcsec -2 . The BCG is indicated by a small 
grey circle. Note that the location of the X-ray centroid is displaced by ~ 4" from the center 
of the gravitationally-lensed arcs. 
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Fig. 2.— RCS0224-0002 Surface Brightness Profiles. Radial surface brightness for the 
0.29-7.0 keV band accumulated in 1" (top panel) and 2" (bottom panel) annuular bins. Solid 
lines trace the best fitting single (3 model fit and dashed lines illustrate the fits obtained by 
fitting this data with double (3 models. The dotted lines represent the best fitting background 
value for each profile. Note the excess emission in the inner 20-30 kpc of the cluster. 
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Fig. 3.— RCS0224-0002 Core Spectra. XSPEC plot showing spectra fit with a single 
temperature spectral model. The spectrum taken from the shorter observation is shown in 
grey. Spectra were extracted from a circle centered on the emission peak with a radius of 
362 hjQ kpc (~ R2500)) an d were grouped to include 20 counts per bin. The fit resulted in 
an ambient cluster temperature of 5.l!o!5 keV. 



Fig. 4. — Hardness Map and Radio Comparisons. Left panel - A hardness map was 
constructed for RCS 0224-0002 by dividing smoothed flux images in the hard (2.0-8.0 keV) 
and soft (0.5-2.0 keV) bands by a smoothed flux image of the total band in the manner Is j~j H . 
Linearly spaced contours were created from the hardness map and are overlaid on a 1.4 GHz 
VLA flux image (gray scale) of RCS0224-0002. Soft emission peaks in a region containing 
both the X-ray centroid (shown as a thick 2.5" radius circle) and a 1.4 GHz radio source. 
Right Panel - HST image of RCS0224-0002 (Figured]) with radio contours overlaid. 1.4 
GHz VLA contours from the flux image shown in the left panel clearly indicate the galaxy 
which is responsible for the radio emission. The position of the X-ray centroid is again 
represented bold 2.5" radius circle. 
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Fig. 5. — Average Tangential Distortion Average tangential distortion is plotted as 
a function of distance from the brightest cluster galaxy, using 462 source galaxies with 
22.5 < F8UW < 25.5 and F606W - F8UW < 1.5 (i.e., bluer than the cluster early types). 
We detect a significant lensing signal. The lower panel shows the signal when the phase 
of the distortion is increased by ir/2. If the signal observed in the upper panel is due to 
gravitational lensing, gx should vanish, as is observed. The solid line is the best fit singular 
isothermal sphere model which yields an Einstein radius of 4.7 ± 0.9 arcseconds. 
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Fig. 6. — Histogram of Galaxy Velocities. Histogram constructed with 39 galaxies, 
indicating a mean redshift of z — 0.778 for RCS0224-0002. The overplotted dashed line 
shows a Gaussian of dispersion 900 km s _1 rest-frame. Note that the two secondary peaks 
at z ~ 0.750 and z ~ 0.805 are visibly separated from the cluster peak and are not spatially 
concentrated, so we do not associate these with the cluster nor with coherent structures. 
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Fig. 7. — X-ray Point Sources. LogN-LogS plot of point sources within R 20 o (1.47 
Mpc) in soft (0.5-2.0 keV) and hard (2-8 ke V) X-ray bands. Th e dashed lines represent 
the X-ray background in these energy bands (IMoretti et al.ll2003l ). The solid line in each 
plot indicates the LogN-LogS distribution of detected sources. Dot dash lines indicate la 
confidence intervals. A slight excess of low flux point sources may be present in the soft 
energy band. 
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Fig. 8. — Point Source Distribution. Histograms of point source density vs. radius were 
constructed from data in soft (0.5-2.0 keV) and hard (2-8 keV) energy bands, and are shown 
here with la error bars. Both energy bands indicate a n excess of point sourc e s with in R,2oo- 
Each of these plots are consistent with the findings of iRuderman fc Ebeling! (120051 ). 



Fig. 9. — VLA Point Source Comparisons. X-ray point source region files produced by 
running WAVDETECT on 0.5-2.0 keV and 2.0-8.0 keV flux images are overlaid on a 1.4 GHz 
VLA image of RCS0224-0002. The size of the point source regions corresponds to the ap- 
parent extent of the sources in the X-ray image. The large circles have a 1.47 h^ Mpc (R200) 
radius. Within that region, at least ten of 26 soft band X-ray point sources are found to be 



Table 1. /3-Model Fits 



Bin Size R cl [h^kpc] fa J i a R c2 [h^kpc] fa J 02 a I B a x 2 /DOF 



1" 


228+1° 


n qqfi+0-069 
u.you_ 058 


fi 7+°- 2 

°-'-0.3 








2 4nq +0 - 007 

z -™ M -0 .007 


223.6/191 




283±i? 


i.i ie±8:8lg 


r O+0.2 
O.O_ 2 


i - ±D -0.07 


U.0/U_ 035 




2 41 1 +0.007 


218.5/188 


2" 


213lJ| 


0.8841°;^ 


7 i +0.5 
' • i -0.4 








2.4i4i8;8i 


153.8/130 




2671^ 


1.03818:8™ 


O.U_ 2 


o 40+O.26 
°-^ y -0.21 


n fi2n+ 0046 

U.OZU_ 038 


fi 3+°- 5 

°-°-0 .4 


2.41618:81 


151.5/127 



Surface brightness I in units of 10 9 photons sec 1 cm 2 arcsec 
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Table 2. Summary of the optical multi-object spectroscopic observations. 



Telescope 


Instrument 


Date 


Spectral Range 


Resolution 


Keck 


LRIS 


11-21-2001 


5800-8000 A 


7.5 A 


VLT 


FORS-2 


01-27-2003 


4500-8500 A 


15 A 






01-28-2003 










01-30-2003 






Gemini 


GMOS-N 


10-21-2003 


4600-9000 A 


16 A 






10-22-2003 










10-23-2003 







-30- 



Table 3. 0.5-2.0 keV Point Sources for R < R 



Number a 


RA 


Dec 


Flux b 




a 


Distance 


1.4 GHz 




[J2000] 


[J2000] 




[10 42 erg s- 1 ] 




[hm Mpc] 




1 


02:24:34.344 


-00:02:26.05 


3.5 ±1.3 


0.61 ±0.23 


3.4 


0.022 




2 


02:24:37.272 


-00:02:49.89 


6.3 ±1.8 


1.11 ±0.31 


5.7 


0.395 


yes 


3 d 


02:24:35.208 


-00:01:28.09 


7.9 ±1.9 


1.39 ± 0.34 


7.7 


0.456 


yes 


4 


02:24:38.208 


-00:03:28.67 


3.2 ± 1.2 


0.56 ±0.22 


3.3 


0.655 


yes 


5(2) d 


02:24:39.024 


-00:01:37.62 


30.2 ±3.7 


5.30 ± 0.65 


27.3 


0.662 




6 


02:24:28.656 


-00:03:18.10 


9.9 ±2.1 


1.75 ± 0.38 


9.5 


0.735 




7(3) 


02:24:33.312 


-00:04:04.98 


5.8 ±1.6 


1.02 ±0.29 


5.9 


0.749 


yes e 


8(4) 


02:24:27.144 


-00:02:25.11 


6.1 ±1.8 


1.07 ±0.32 


5.4 


0.794 


yes 


9 


02:24:38.112 


-00:00:56.40 


3.2 ± 1.3 


0.56 ±0.22 


3.1 


0.813 




10 (5) d 


02:24:28.728 


-00:03:42.58 


32.8 ±3.8 


5.76 ± 0.67 


28.2 


0.841 




11 


02:24:27.288 


-00:01:32.81 


5.7 ± 1.7 


1.00 ±0.29 


5.3 


0.876 


yes 


12 (8) d 


02:24:40.368 


-00:03:40.72 


61.3 ±5.2 


10.78 ±0.92 


56.5 


0.898 


yes 


13 (7) 


02:24:42.048 


-00:02:43.56 


18.5 ±2.9 


3.25 ±0.51 


17.9 


0.901 




14 (9) 


02:24:33.984 


-00:04:26.01 


24.8 ± 3.3 


4.36 ± 0.58 


23.7 


0.902 




15 (10) 


02:24:30.600 


-00:00:38.78 


1429.4 ± 25.3 


251.23 ±4.44 


916.0 


0.906 


yes 


16 (11) 


02:24:26.112 


-00:01:54.61 


8.0 ±1.9 


1.40 ± 0.34 


7.6 


0.941 


? 


17 


02:24:39.648 


-00:00:46.28 


5.3 ±1.6 


0.93 ± 0.28 


4.8 


0.977 


yes 


18 (12) 


02:24:28.128 


-00:04:03.81 


5.3 ±1.6 


0.93 ± 0.28 


5.4 


1.002 


yes 


19 (13) 


02:24:40.176 


-00:04:09.33 


39.5 ±4.2 


6.95 ± 0.74 


36.9 


1.031 




20 (14) 


02:24:42.384 


-00:03:32.32 


4.5 ± 1.5 


0.79 ±0.26 


4.5 


1.055 




21 


02:24:25.200 


-00:03:11.21 


5.4 ± 1.6 


0.94 ± 0.28 


5.1 


1.069 




22 


02:24:29.544 


-00:00:20.37 


29.2 ±3.7 


5.13 ±0.64 


25.2 


1.084 




23 


02:24:43.032 


-00:03:26.71 


7.6 ± 1.9 


1.34 ±0.33 


7.7 


1.101 




24 (16) 


02:24:38.784 


-00:00:05.81 


8.0 ±2.0 


1.41 ±0.36 


6.9 


1.183 




25 (18) 


02:24:38.160 


+00:00:09.60 


5.6 ±1.6 


0.98 ±0.29 


5.4 


1.260 




26 d 


02:24:21.456 


-00:02:08.01 


3.6 ± 1.4 


0.64 ± 0.24 


3.4 


1.442 





a 0.5-2.0 keV (2.0-8.0 keV) 

b Unabsorbed 0.5-2.0 keV flux in units of 10~16 erg cm -2 s _1 

c Unabsorbed bolometric X-ray luminosity obtained using 0.5-2.0 keV count rates and PIMMS, assum- 



31 



ing z=0.778 and a powerlaw index of 1.6 

d Within 1" of a possible cluster member, using photometrically determined redshifts 



-SCUBA submillimeter source (jWebb et al.ll2005l ) 
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Table 4. 2.0-8.0 keV Point Sources for R < R 



200 



Number a 


RA 
[J2000] 


Dec 
[J2000] 


Flux b 


L x c 

[10 42 erg s" 1 ] 


a 


Distance 
[h-J Mpc] 


1.4 GHz 


1 


02:24:36.408 


-00:02:03.15 


19.4 ± 6.9 


1.96 ± 0.70 


3.5 


0.309 


yes 


2 (5) d 


02:24:39.000 


-00:01:37.52 


34.3 ± 8.8 


3.46 ± 0.89 


6.4 


0.660 




3 (7) 


02:24:33.336 


-00:04:04.34 


36.9 ± 8.9 


3.72 ± 0.90 


7.3 


0.744 


yes e 


4 (8) 


02:24:27.120 


-00:02:25.40 


107.5 ± 16.0 


10.85 ± 1.62 


16.1 


0.796 


yes 


5 (10) a 


02:24:28.728 


-00:03:42.14 


28.0 ± 8.0 


2.83 ± 0.81 


5.3 


0.838 




6 d 


02:24:40.944 


-00:01:24.94 


19.5 ± 6.8 


1.97 ± 0.69 


3.8 


0.897 




7(13) 


02:24:42.024 


-00:02:43.26 


36.1 ±9.0 


3.64 ±0.91 


6.7 


0.898 




8 (12) d 


02:24:40.368 


-00:03:40.75 


119.8 ± 15.9 


12.10 ± 1.61 


19.9 


0.899 


yes 


9 (14) 


02:24:33.984 


-00:04:25.94 


37.0 ±8.9 


3.73 ± 0.90 


7.5 


0.901 




10 (15) 


02:24:30.600 


-00:00:38.86 


862.3 ±42.2 


87.05 ± 4.26 


119.0 


0.906 


yes 


11 (16) 


02:24:26.112 


-00:01:54.13 


16.1 ±6.1 


1.63 ±0.62 


3.4 


0.942 




12 (18) 


02:24:28.104 


-00:04:04.18 


20.3 ±6.9 


2.05 ±0.70 


3.9 


1.007 


yes 


13 (19) 


02:24:40.176 


-00:04:09.44 


74.1 ± 12.6 


7.48 ± 1.27 


13.6 


1.032 




14 (20) d 


02:24:42.384 


-00:03:31.73 


37.4 ±9.2 


3.77 ±0.93 


6.6 


1.053 




15 


02:24:28.824 


-00:04:38.94 


20.0 ±7.3 


2.02 ±0.73 


3.4 


1.167 




16 (24) 


02:24:38.856 


-00:00:06.35 


27.5 ±8.4 


2.78 ±0.85 


4.6 


1.182 




17 d 


02:24:37.344 


-00:04:57.77 


124.6 ± 16.3 


12.58 ± 1.64 


20.4 


1.196 




18 (25) 


02:24:38.160 


+00:00:09.53 


74.0 ± 12.9 


7.47 ± 1.30 


12.6 


1.260 





a 2. 0-8.0 keV (0.5-2.0 keV) 

b Unabsorbed 2.0-8.0 keV flux in units of 10~16 erg cm -2 s _1 

c Unabsorbed bolometric X-ray luminosity obtained using 2.0-8.0 keV count rates and PIMMS, as- 
suming z=0.778 and a powerlaw index of 1.6 

d Within 1" of a possible cluster member, using photometrically determined redshifts 



"SCUBA submillimeter source (jWebb et al.ll2005l ) 



